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F&R Docket No.: 07508-046002 

A Gain Factor and Position Determination System 

TECHNICAL FIELD 

This invention relates to systems that use magnetic fields to determine an object's 
location and orientation, and system gain factor. 

BACKGROUND 

As is known in the art, systems may use magnetic field measurements to indirectly 
determine the location and orientation of an object. These systems are useful, for example, in 
the medical field, because they are able to accurately locate an object in a patient's body with 
only minimal intrusion into the body. The intrusion involves placing a small probe near the 
object to be located. The three-dimensional location and orientation of the probe is then 
determined from the effect that the probe's location and orientation have on magnetic field 
measurements. 

The probe may be either a source or a sensor of a magnetic field. If the probe is a 
source, sensors exterior to the body measure the field produced by the probe. If the probe is 
a sensor, magnetic sources exterior to the body produce the fields being measured. 

Determining a probe's location and orientation from magnetic field measurements is 
not straight forward because the measured magnetic fields are nonlinear functions of the 
location and orientation. To determine the probe's location and orientation from the 
measured magnetic field values, the probe's location and orientation are first presumed or 
"guessed" to be at a predicted location and orientation. An iterative process is used to 
compare values of the magnetic field at the guessed probe location and orientation with the 
measured field values. If the magnetic field values at a guessed location and orientation are 
close to the measured values, the guessed location and orientation are presumed to accurately 
represent the actual location and orientation of the probe. 

The iterative process uses a physical model for the probe's environment. The 
physical model specifies the location and orientation of each field source. From the specified 
locations and orientations, laws of electrodynamics determine the field values. 



1 



F&R Docket No.: 07508-046002 

As the probe and its positioning system are physical systems, they are susceptible to 
various external influences (e.g., stray magnetic fields, field absorbing materials being 
positioned proximate the field generators and/or sensors, etc.) that affect the gain of the 
system. Additionally, these physical devices have various engineering tolerances (e.g., cable 
5 resistance, probe gain, input impedance, etc.) that also affect overall system gain. 

Accordingly, each time a component of the system is replaced, the system must be manually 
recalibrated. 

SUMMARY 

10 According to an aspect of this invention, a system for determining the position, 

orientation and system gain factor of a probe includes magnetic field sources and at least one 
magnetic field sensor, such that a combination of a magnetic field sensor and a magnetic 
field source generates a unique measured magnetic field value. The system also includes a 
probe whose position and orientation affect the unique measured magnetic field values. A 

15 processor, coupled to receive these unique measured magnetic field values, iteratively 

processes measured magnetic field values to determine a system gain factor indicative of the 
gain of the probe and location factors indicative of the position and orientation of the probe. 
The number of unique measured magnetic field values generated must be at least equal to the 
sum of the number of factors calculated. 

20 One or more of the following features may also be included. The iterative process is 

configured to determine a function of the differences between the measured magnetic field 
values and a plurality of predicted magnetic field values. The processor includes a calculated 
location process for calculating the predicted magnetic field values, in that the calculated 
location process guesses an initial gain, position and orientation for the probe, and then 

25 calculates the predicted magnetic field values based on a physical model and the initial gain, 
position and orientation. The initial position and orientation may be a predetermined or 
randomly selected fixed point. 

The processor includes an optimization function for determining an extremum 
indicative of the differences between the measured magnetic field values and the predicted 

30 magnetic field values. The optimization function is a least squares sum function. The 
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processor includes a repositioning process for adjusting the initial gain, position and 
orientation of the probe in response to the extremum being in a predefined range of 
unacceptable values, which is indicative of an unacceptable level of difference between the 
measured magnetic field values and the plurality of predicted magnetic field values. The 

5 location factors may include spatial, spherical, and/or rotational coordinates. 

According to a further aspect of this invention, a method for determining the position, 
orientation and system gain factor of a three-dimensional object includes positioning a 
plurality of magnetic field sources proximate the three-dimensional object and positioning at 
least one magnetic field sensor in a fixed spatial relationship with that three-dimensional 

10 object. A combination of a magnetic field sensor and a magnetic field source generates a 
unique measured magnetic field value. Further, the position and orientation of the three- 
dimensional object affects these unique measured magnetic field values. The method 
determines a system gain factor indicative of the gain of the three-dimensional object and a 
plurality of location factors indicative of the position and orientation of the three-dimensional 

15 object. The number of unique measured magnetic field values generated must be at least 
equal to the sum of the number of factors calculated. 

One or more of the following features may also be included. The step of determining 
a system gain factor and a plurality of location factors includes determining a function of the 
differences between the measured magnetic field values and a plurality of predicted magnetic 

20 field values. The step of determining a system gain factor and a plurality of location factors 
includes guessing an initial gain, position and orientation for the three-dimensional object 
and calculating the predicted magnetic field values based on a physical model and the initial 
gain, position and orientation. The step of determining a system gain factor and a plurality of 
location factors includes determining an extremum indicative of the differences between the 

25 measured magnetic field values and the predicted magnetic field values. The step of 

determining a system gain factor and a plurality of location factors includes adjusting the 
initial gain, position, and orientation of the three-dimensional object in response to the 
extremum being in a predefined range of unacceptable values, which is indicative of an 
unacceptable level of difference between the measured magnetic field values and the 

30 plurality of predicted magnetic field values. 
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While the system and method described above includes a plurality of magnetic field 
sources and at least one magnetic field sensor, the system can also include a plurality of 
magnetic field sensors and at least one magnetic field source. Further, while the system and 
method described above were said to include a probe, that probe can actually be any three- 
5 dimensional object, such as a hollow tube (e.g., a biopsy needle). 

The advantages of the above aspects of the invention are numerous. As this system 
automatically calculates the system gain factor for the probe or three-dimensional object 
being utilized, device gain calibration is automated. Therefore, the need for tedious and 
time-consuming manual gain recalibration is eliminated. Additionally, as this calibration 
10 process is automated, the system components (e.g., probes, sensors, leads, etc.) can be 

quickly and easily swapped and the system reconfigured without the need for manual gain 
recalibration. This, in turn, streamlines and simplifies the reconfiguration process. 
Additionally, as the system automatically and continuously determines the system gain 
factor, the system can be physically moved without fear of external influences mandating 
1 5 manual gain recalibration. 

The details of one or more embodiments of the invention are set forth in the accompa- 
nying drawings and the description below. Other features, objects, and advantages of the 
invention will be apparent from the description and drawings, and from the claims. 

DESCRIPTION OF DRAWINGS 

20 FIG. 1 is a schematic representation of a system that uses magnetic measurements to 

find a probe's system gain factor, location, and orientation; 

FIG. 2 is a perspective view of the system of FIG. 1; 

FIG. 3 is a side view of an alternate system that uses magnetic measurements to find a 
probe's gain, location, and orientation; 

25 FIG. 4 is a flow chart for a process that iteratively guesses the gain, location, and 

orientation of a probe from measured magnetic field values; 

FIG. 5 is a flow chart for a process that determines whether the best guess of the gain, 
location, and orientation of a probe is reliable; 
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FIG. 6 is a flow chart for a process that evaluates an error or optimization function at 
global and local extrema; 

FIG. 7 is a flow chart for a process that determines extrema values of the error or 
optimization function that correspond to the presence of distortions; and 

5 FIG. 8 shows a computer that finds a probe's gain, location, and orientation from 

measurements of magnetic fields. 

Like reference symbols in the various drawings indicate like elements. 

DETAILED DESCRIPTION 

Referring to FIGs. 1 and 2, there is shown a schematic representation of a system 10 

10 that uses measurements of magnetic fields to find the location, orientation and system gain 
factor of a moveable probe 12. The moveable probe 12 is located inside a volume 14 (e.g., a 
body of a medical patient). System 10 also includes a plurality of field sources 16i_ n (e.g., 
small induction coils) that are located outside of volume 14. Field sources 16i_ n are driven by 
a field-generating module 22 (e.g., an alternating or direct current source). Please note that 

15 while seven field sources are shown, this is for illustrative purposes only and is not intended 
to be a limitation of the invention. 

FIG. 2 is a side view of system 10 that shows the three-dimensional orientation and 
location of field sources 16i_ n and probe 12. Each field source 16i- n is either: a surface- 
mount field source 17 located on a surface 18 of a regular tetrahedron 24; or an edge-mount 

20 field source 19 located on an edge 20 of a regular tetrahedron 24. The field sources 16i- n are 
typically oriented so that their internal magnetic moments "m" are either perpendicular to the 
surfaces 18 of tetrahedron 24 or parallel to the edges 20 of tetrahedron 24. Probe 12 is 
located outside of tetrahedron 24 and has an orientation defined by a normal direction "n" to 
a sensor 13 that measures local magnetic field values. Sensor 13 may include a simple coil, 

25 several coils, a Hall sensor, or a flux gate sensor and may measure a magnetic field flux or 
magnetic field differential. Note that while system 10 is shown to include a single 
tetrahedron 24 having a maximum of four surface-mount field sources 17 and six edge- 
mount field sources 19, this is for illustrative purposes only and is not intended to be a 
limitation of the invention, as the positioning of these field sources 16i_ n can be adjusted and 

30 reconfigured as known by those skilled in the art. For example, a pair of tetrahedrons 24 
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utilizing only edge-mount field sources 19 can be utilized. This configuration allows for a 

maximum of twelve distinct field sources. 

Each source 16i_ n may also include single or multiple field coils. For a source 16i_ n 

with a single coil, the magnetic field in the volume 14 can be approximated by a dipole field. 
5 For a source 16i. n with several coils, the magnetic field in the volume 14 may be a higher 

multipole field. In one embodiment, each source 16i_ n uses two identical coils whose 

normal vectors are anti-parallel. Such a source produces a magnetic quadrapole field in 

volume 14. A quadrapole field has larger spatial variations than a dipole field and thus, may 

be more convenient for locating probe 12. 
10 Other embodiments may use field sources 16i_ n that primarily generate higher 

magnetic multipole fields. Further, other embodiments may use a different numbers of field 

sources 16i- n or locate and orient the field sources 16i_ n differently. 

Referring to FIGs. 1 and 2, magnetic fields from sources 16i_ n induce electromotive 

forces (EMF f s) in internal sensor coil 13 of probe 12. An electronics module 26 connected to 
15 probe 12 measures the EMF's. The measured EMF's represent the measured local values of 

magnetic fields at the probe's location and orientation in space 14. Electronics module 26 

also can identify the individual source 16i_ n that produces each measured EMF. 

In one embodiment, measurement timing information is used to identify the source 

16i- n producing a measured field. In this embodiment, field-generating module 22 
20 temporally multiplexes power to the different field sources 16i_ n and relays timing 

information to the electronics module 26. 

In another embodiment, the field generating module 22 drives each field source 16i_ n 

at a different frequency. To identify the particular source 16i_ n of a measured field, 

electronics module 26 or computer 28 decomposes measured EMF's from the probe's coil 13 
25 into frequency components. These frequency components of the measured fields are then 

matched to individual field sources 16i_ n . 

In either embodiment, electronics module 26 outputs several measured magnetic 

fields Bi. n measured corresponding to the individual field measurements performed. It is 

important to note that there are several hardware configurations that can generate these field 
30 measurements. For example, if eight distinct field measurements are desired, system 10 can 

be configured so that eight field sources 16i_ n are utilized in conjunction with a probe 12 
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incorporating one sensor coil 13. In this configuration sensor coil 13 would measure the field 
generated by each of the eight field sources 16i_ n , resulting in eight distinct field 
measurements. 

Alternatively, a probe 12 could be constructed that incorporates multiple sensor coils 
5 13. For example, if a probe 12 was constructed that includes two sensor coils 13 

(configuration not shown), this configuration would double the number of individual field 
measurements achievable by that probe (using the equivalent number of field sources 16i- n ). 
Specifically, if probe 12 includes two sensor coils 13, each coil could independently measure 
the strength of the magnetic field generated by a single field source. Therefore, if (as stated 

10 above) eight distinct field measurements are desired and a probe 12 is utilized that includes 
two sensor coils 13, only four field sources 16i_ n would be required, as each coil 13 could 
independently measure the field generated by each of the four field sources, thus resulting in 
eight distinct field measurements. 

These measured magnetic field values Bi_ n measured depend on the system gain factor of 

15 probe 12 and the three-dimensional location and the orientation of the probe's coil 13. Please 
note that the lower case "n" is used to denote the number of measured magnetic fields 
(except where it denotes normal vectors in the various figures). 

As stated above, the number of field sources 16i_ n and sensor coils 13 incorporated 
into and utilized by system 10 varies depending on the number of factors to be determined by 

20 the system. Specifically, system 10 will determine the system gain factor of probe 12. 

Additionally, system 10 will also determine the position and orientation of probe 12. As the 
position and orientation of probe 12 is described by specifying up to six degrees of freedom 
(i.e., x-axis position, y-axis position, z-axis position, roll, pitch, and yaw) for probe 12, the 
maximum number of position factors that system 10 will calculate is six. Therefore, the 

25 maximum number of factors (both positional and gain) that system 10 calculates is seven. 
Thus, the number of distinct field measurements required to determine these factors is one 
greater than the number of factors being determined. Accordingly, if system 10 is 
determining the system gain factor plus six positional factors (i.e., degrees of freedom), a 
total of seven calculated factors need to be determined. This would require eight distinct 

30 field measurements. As stated above, this can be achieved utilizing a probe 12 incorporating 
a single sensor coil 13 and eight field sources 16i_ n . Alternatively, this could be achieved by 
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utilizing a probe 12 incorporating eight sensor coils 13 (if physically possible) and only one 
field source 16i_ n . Accordingly, the important issue is the number of distinct field 
measurements made and not the specific number of field source 16i_ n or sensor coils 13 
employed. 

5 Similarly, if system 10 is determining the system gain factor plus five positional 

factors (i.e., five degrees of freedom), a total of six calculated factors need to be determined. 
Again, as above, this can be accomplished utilizing a variety of configurations. 

As will be described below, the use of a number of distinct field measurements that is 
greater than the number of factors being calculated allows for the determination of the quality 

10 of the guesses used during an iterative process. If this quality determination is not required 
(or desired), the specific number of distinct field measurement required would be equal to the 
number of factors being calculated. 

Further, when employing multiple field sensors, it is important to realize that there is 
a gain factor associated with each sensor. Accordingly, for the example stated above 

15 concerning the determination of six degrees of freedom, since at least two non-coaxial sensor 
coils are required to determine six degrees of freedom, there are at least two system gain 
factors (one for each sensor coil). Therefore, when calculating the six degrees of freedom 
and system gain factor of a probe utilizing two non-coaxial sensor coils, the system gain 
factor actually consists of two components, namely a gain factor for each sensor within the 

20 probe. Accordingly, the total number of calculated factors is eight (six positional, the gain 
factor for the first field sensor, and the gain factor for the second field sensor) and, therefore, 
the number of distinct field measurements required is nine (if the quality determination is 
required or desired). 

Electronics module 26 sends the measured field values to computer 28. Computer 28 
25 uses the measured magnetic field values to determine the probe's system gain factor and 
location / orientation by comparing the measured magnetic field values to magnetic field 
values from a physical model, as will be described below in greater detail. 

The physical model is a set of physical equations that determine values of magnetic 
fluxes measured by probe 12 as a function of several parameters. The parameters include: 
30 the position, orientation, and magnetic moments of field sources 16i_ n ; the location, 

orientation, and sensitivity of probe 12; and characteristics of electronics module 26. A 
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vector (x, y, z) and a pair of angles (<p, 0) specify the three-dimensional location and 
orientation of sensor coil 13 of probe 12. If the probe 12 has multiple non-colinear coils, the 
parameters may include an additional angular parameter (\|/) that defines a rotational aspect 
of the probe 12. Note that this factor (i.e., the sixth degree of freedom) can only be 
5 calculated by utilizing a probe having a second coil on a different axis (non-coaxial), as 

multiple coils about the same axis would not allow system 10 to sense probe rotation about 
that axis. 

As is known by those skilled in the art, the physical model may describe each source 
as a magnetic multipole so that the fields measured by sensor coil 13 are the associated 

10 multipole fields (e.g., dipole or quadrapole). The multipole field values depend on the 

system gain and the location, orientation, and magnetic moment "m" of each source 16i_ n . 
The measured values of the magnetic flux depend on the location, size, orientation and gain 
of sensor coil 13 with respect to field sources 16i. n . 

The physical model is also based on underlying assumptions about the environment 

15 near the volume 14. For example, the model assumes preselected values for the location and 
orientation of each field source 16i_ n and the absence of other sources or field distorting 
objects. The presence of field distorting objects 30, 32 (e.g., conductors or new field 
sources) may invalidate the model's predictions for field values. However, the negative 
effects of constant background fields are eliminated, as sensing coil 13 only measures time 

20 varying magnetic fields. Alternatively, if static field measurements are desired, a flux gate 
sensor or hall effect sensor (as is known in the art) can be utilized, as they allow for 
measurement of static (or constant) magnetic fields. 

FIG 3 shows another system 40 that uses magnetic measurements to find the system 
gain factor and location / orientation of a moveable probe 42. In system 40, the roles of the 

25 field sensors and sources are interchanged. The movable probe 42, which is located inside an 
observation volume 44, is a source of a magnetic field and external field sensors 46i_ n 
measure the magnetic field produced by probe 42. Probe 42 connects to a field-generating 
module 51 (e.g., a voltage source). As above, these field sensors 46i. n , which are located on 
either the edges or surfaces of a tetrahedron 41, measure either magnetic fields or magnetic 

30 field gradients through induced EMF's. Each sensor 46i_ n may have one or more magnetic 
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field sensors oriented to measure fields in various directions. Each sensor 46i_ n has an 
orientation "p" fixed by the orientation of its one or more internal magnetic field sensors. 

Please note that, as stated above, the location and orientation of each sensor, as well 
as the number of sensors 46i- n utilized, may vary depending on the specific application and 
5 the number of factors (i.e., position / orientation and system gain) being determined. 

An electronics module 52 monitors EMFs from the various field sensors 46i- n to 
determine the individual magnetic field values. These measured magnetic field values are 
then sent to a computer 54, which calculates the system gain factor and location / orientation 
of the probe 42 from these measured magnetic field values. 

10 Referring to FIGs. 1-3, both systems 10, 40 measure a set of magnetic fluxes to obtain 

a set of measured magnetic field values Bi_ n measured ? such that "n" is greater than the number 
of factors (i.e., position and system gain) being calculated. 

This set of measured field values Bi. n measured , which are obtained from measured 
magnetic fluxes, also has a non-linear dependence on the three-dimensional location / 

15 orientation of the probe and a linear dependence on the system gain factor. As stated above, 
the probe's location and orientation are defined by a vector (x, y, z) and at least a pair of 
azimuthal and polar angles (cp, 0), respectively. Further, the system gain factor of probe 12 is 
defined by a gain coefficient (g). By using a physical model for the "measured" field 
dependencies, the systems 10, 40 can iteratively determine the probe's gain factor, location, 

20 and orientation from the associated set of measured field values Bi_ n measured . 

The physical model describes a preselected magnetic environment in the region of the 
field sensor[s] (i.e., the volumes 14 and 44 shown in FIGs. 1 and 3 respectively). The 
preselected magnetic environment may or may not include contributions from nearby 
conducting objects (i.e., objects 30 and 32). If the preselected environment is different from 

25 the actual environment, the model may predict incorrect magnetic field values; otherwise the 
predictions are correct. The actual environment may be different due to the presence of field 
distorting objects 30, 32. Field distorting objects 30, 32 include conducting objects that 
support eddy currents (e.g., a pair of surgical scissors, ferromagnetic materials, and active 
sources of magnetic fields). The presence of such objects can invalidate magnetic 

30 determinations of the probe's gain, location, and orientation. 
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The iterative process may also yield incorrect probe gain, location, and orientation 
because of hardware or software failures in modules 22, 26, 51, 52 or computers 28, 54. 

The presence of distorting conditions may not be evident to users of systems 10, 40. 
The users may be unskilled with respect to physical basis for the magnetic location systems 
5 (e.g., the users may be medical doctors). To avoid errors by unskilled users, each system 10, 
40 detects and warns users about the presence of potentially measurement distorting 
conditions (e.g., by flashing messages on a video monitor or through audio alert signals). 

Please note that while probe 12, 42 is stated as being a probe, this is only intended to 
ease in the explanation of the subject invention and is for illustrative purposes only. Probe 12 

10 can be various other devices known in the art (e.g., a catheter, an endoscope, biopsy needles, 
body-mounted position sensors, etc.). 

FIG. 4 shows a flow chart for an iterative process 60 that uses the measured magnetic 
field values Bi_ n measured to determine the gain, location, and orientation of either probe 12 of 
FIGs. 1-2 or probe 42 of FIG. 3. 

15 Process 60 receives 62 an initial guess for the probe's system gain factor, location 

and orientation. This initial guess can be predefined in the process, user definable, or 
randomly selected. The initial guess is a preselected point of the (x, y, z, (p, 6) parameter 
space that defines the location and orientation of the probe and a preselected system gain 
factor (g). The initial guess is the first accepted guess for the probe's gain, location, and 

20 orientation. From the last accepted guess, the process makes a new guess 64 for the probe's 
location, orientation, and gain. Naturally, while the initial preselected point (x, y, z, cp, 0) is 
shown to include five degrees of freedom, this is for illustrative purposes only and is not 
intended to be a limitation of the invention. Specifically, if greater or fewer degrees of 
freedom are desired, the preselected point (x, y, z, (p, 0) will have greater or fewer variables 

25 respectively. 

Each new guess of the probe's gain, location, and orientation is found from the last 
accepted guess by any of a variety of procedures including the Levenberg-Marquardt 
process, a log-likelihood function, a neural network, simulated annealing, a genetic 
algorithm, a simplex process, or other processes known to those skilled in the art 
30 The Levenberg-Marquardt process is an iterative process used to find a best match 

between a set of measurements and a set of values obtained from preselected nonlinear model 
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equations. The Levenberg-Marquardt process is described in: Numerical Recipes in C: the 
Art of Scientific Computing , by W.H. Press et al., Cambridge University Press 1992 and is 
incorporated herein by reference. 

In Levenberg-Marquardt process, the model equations are the set of physical 
5 equations B,. n predicted (x, y, z, <p, 0, g), which define the magnetic field values predicted in 
terms of the probe's location, orientation, and gain (x, y, z, cp, 0, g). The model equations 
come from the physical laws of electrodynamics. In one embodiment, the model equations 
may describe the magnetic field of each field source as a magnetic dipole or as a magnetic 
quadrupole. 

10 The Levenberg-Marquardt procedure iteratively tries to find a best match between the 

measured magnetic field values Bi_ n measured and the magnetic field values B,_ n predicted 
predicted from the physical model equations. The N th accepted guess for a match is 
associated with the probe gain, location, and orientation coordinates (x N , y N , zn, <Pn, 0n> gN). 
From these coordinates and the physical model equations, the Levenberg-Marquardt process 

15 generates an (N+l) th guess that specifies the associated gain, location, and orientation 
coordinates (x N +i, Yn+u zn+i; 9n+i> 0n+i> gN+i) of the probe. The Levenberg-Marquardt 
equation for the (N+l) th guess utilizes the values of the fields B Kn predicted and derivatives of 
the fields Bi_ n P redicted evaluated at the value of the N th accepted guess for a match. The 
Levenberg-Marquardt process provides new guesses that rapidly produce the best match 

20 between measured field values Bi_ n measured and predicted field values Bi_ n predicted which are 
obtained from the non-linear model equations. Please note that the lower case "n" is used to 
denote the number of measured magnetic fields (e.g., B 1 . n measured ) and the upper case "N" is 
used to denote the specific guess processed by the Levenberg-Marquardt procedure (e.g., 
(N+l th )). 

25 The process 60 evaluates the quality of each new guess for the probe's gain, location, 

and orientation. To determine the quality of this guess, process 60 calculates 66 new 
magnetic field values which correspond to the new guess for the probe gain, location and 
orientation (i.e., Bl . n predicted(new) = B 1 . n predicted (x N+1 , y N+l , z N+1 ; <p N+1 , 0 N+1 , g N+1 ) for the 
(N+l) th guess. Using both the calculated (i.e., B Un predicted ) and the measured (i.e., Bi. n measured ) 

30 magnetic field values, process 60 evaluates 68 an optimization function. The optimization 

function is also referred to as an error function, because the optimization function is sensitive 
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to differences between measured (i.e., Bi_ n measured ) and calculated (i.e., B,. n predicted ) field 
values. A global extremum of the optimization function defines the "best" guess for the 
location and orientation of the probe. 

One embodiment uses a least-squares sum (i.e., x 2 ) as the optimization function. In 
5 general, these functions have both minima and maxima extrema. However, the extrema of 
the least-squares sum that we are interested in are the minima. The value of the least-squares 
sum x 2 (N) for the magnetic field associated with the N th guess has the form: 

X 2 (N) = Zi-n [Bi- n measured - B 1 ./ redicted (x N , y N , z N , cp N , 0 N , g N )] W. 

The sum goes over the "1-n" members in the set of measured field values (i.e., Bi_ n 

10 measured^ ^ ^ 0 b ta i ne( j f or a s i ng i e p ro be location and orientation. The term c\. n represents 
uncertainties associated with the measurement of Bi_ n measured . 

Process 60 determines 70 whether the value of the optimization function for the new 
guess is closer to a value at an extremum than the value for the last accepted guess. For the 
least-squares sum, since the extrema are minima, the new value is closer to a minimum if 

15 x (N+l) < X (N). If the value of the optimization function for the new guess is closer to an 
extremum, process 60 accepts 72 the new guess for the gain, location, and orientation of the 
probe. If the value of the optimization function for the new guess is farther from an 
extremum (e.g., x 2 (N+l) > % 2 (N))> process 60 rejects 74 the new guess. After accepting 72 
the new guess, process 60 determines 75 if the difference between the optimization function 

20 for the old and new guesses is less than an acceptable threshold (e.g., 0.01). If it is less than 
this threshold, the probe's gain, location, and orientation are reported 77. If the new guess is 
rejected 74 or if the difference between the optimization function for the old and new guesses 
is not less than an acceptable threshold, process 60 increments 76 a counter for the number of 
iterations performed. Process 60 then checks 79 to see if the incrementation counter is 

25 greater than the maximum number of iterations specified. If it is, process 60 then reports 77 
the probe's gain, location, and orientation. If the number of iterations performed is not 
greater than the maximum number of iterations specified, process 60 loops back 78 to find a 
new and better guess of the probe's gain, location, and orientation. 

Process 60 outputs the last accepted guess for the probe's gain, location, and 

30 orientation and a count specifying the number of iterations performed. In some 
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embodiments, process 60 performs a preselected number of iterative loops 78 to find a better 
guess before reporting the accepted guess for the probe's gain, location, and orientation. This 
produces reported guesses that are closer to values associated with extremum of the 
optimization function. 

5 FIG. 5 is a flow chart for a process 80 that uses magnetic field measurements to 

determine a probe's gain, location, and orientation. Process 80 provides an initial guess 82 
for the probe's gain, location, and orientation. The initial guess of the probe's gain, location, 
and orientation (x, y, z, (p, 9, g) may include either a preselected fixed point or a randomly 
selected point in the (x, y, z, (p, 9) space. For the selected initial guess, process 80 performs 

10 84 iterative process 60 (of FIG. 4) to obtain a better guess for the probe's gain, location, and 
orientation. The better guess is based on the selected initial guess and the measured magnetic 
field values. Process 80 may perform several iterative loops of process 60 to obtain a better 
guess, which is closer to an extremum of the optimization function used in process 60. The 
optimization process 60 provides a value for the optimization function (e.g., the % 2 function) 

15 and the loop count (which indicates the number of iterations performed to obtain the better 
guess). 

The value of the optimization function provides data indicative of the reliability of the 
better guess of the probe's gain, location, and orientation. Random measurement errors cause 
the value of the optimization function to fall on a probabilistic distribution function whose 

20 form is independent of the physical model for the magnetic field measurements. For the 
least-squares sum, the probabilistic distribution function is known as the x 2 -distribution. 
Systematic measurement errors also affect the value of the optimization function. 

The extrema of the optimization function can be maxima or minima and can fall into 
several classes. An extremum may be local or global extremum. Local and global extrema 

25 are distinguished by the associated values of the optimization function. For the least-squares 
sum, the value of the optimization function at a global minimum is smaller than the value of 
the function at a local minimum. Thus, global and local minima are associated with 
respective low and high values of the least-squares sum. 

An extremum may also correspond to a situation in which measurements of magnetic 

30 fields are either distorted or undistorted. For the least squares sum, values of the 

optimization function at global minima for undistorted measurements are smaller than values 
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of the optimization function at global minima for distorted measurements. Also, values of 
the optimization function at global minima for distorted measurements are smaller than 
values of the optimization function at local minima. 

Thus, values of the optimization function at extrema carry information about the 
5 quality of the estimates of a probe's gain, location, and orientation obtained through process 
60. The information enables determining whether random or systematic errors are present. 
Values of the least- squares sum at extrema are typically ordered. The lowest values 
correspond to global minima in the absence of distortions to field measurements. 
Intermediate values correspond to global minima in the presence of distortions to field 

10 measurements. The highest values correspond to false or local minima for which estimates 
of the probe's gain, location, and position are not reliable. 

Distortions may occur during generation of magnetic fields, measurement of 
magnetic fields, acquisition of field measurements, or processing of field measurements. The 
distortions to the generation of magnetic fields may occur due to a failure of a field source 

15 16i_ n , probe 42, or field-generating module 22 or 51 . The distortion of magnetic field 

measurements may result from the presence of conductors or ferromagnetic materials in the 
monitored space, which distort the time-varying magnetic fields. Distortion during the 
acquisition or processing of magnetic field measurements may result from hardware or 
software failures (e.g., in electronics modules 26, 52 or computers 28, 54). 

20 Prior to performing process 80, a calibration is performed to classify extrema of the 

optimization function. The calibration classifies extrema of the optimization function into 
three or more sets (Si). One set (Sg-nd) corresponds to real global extrema of the 
optimization function for which field measurements and field measurement processing is not 
distorted. Another set (Sl) corresponds to false or local extrema of the optimization function. 

25 A third set (Sg-d) corresponds to real global extrema of the optimization function for which 
either field measurements or field measurement processing is distorted. 

New sets may be formed from the sets (Sg-nd, Sl, and Sg-d) by set intersection and 
union operations. One set (Snd) includes values of the optimization function that are only 
associated with extrema of the function for which field measurements and measurement 

30 processing are not distorted. This set (S N d) is defined as Sg-nd - (Sl u Sg-d). Another set 

(S D ) includes values of the optimization function that are only associated with global extrema 
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of the function for which field measurements or field measurement processing is distorted. 
This set (So) is defined as Sg-d - (Sl u Sg-nd). Another set (Slo) includes values of the 
optimization function that are only associated with false or local extrema of the function. 
This set (Slo) is defined as Sl - (Sg-nd u Sg-d). Finally, a set (Snd-d) includes values of the 
5 optimization function that are only associated with global extrema of the function. For 

values in Snd-d, the field measurements and field measurement processing may or may not be 
distorted. Snd-d is defined as Sg-nd u Sg-d- In various embodiments, some of the above- 
described sets (Si) may be empty. 

Referring again to FIG. 5, process 80 uses the calibrated classification of values of the 

10 optimization function to classify extrema found with iterative process 60. Process 80 

determines 86 whether the value of the optimization function at each extremum corresponds 
only to a global minimum without distortion (i.e., whether the value belongs to Sg-nd). If the 
value belongs to set Sg-nd, process 80 registers 88 the associated "better" guess for the 
probe's gain, location, and orientation as the probe's gain, location, and orientation. For 

15 example, the "better" guess for the location / orientation coordinates and the gain (x, y, z, (p, 
0, g) may be displayed on a computer screen for a user to view as the final estimate of the 
gain, location, and orientation of the probe. Further, process 80 determines 87 if the system 
gain factor is within an acceptable range. If so, process 80 stops. If not, process 80 provides 
93 a warning to the user that the system gain factor is out of range and, therefore, the probe's 

20 operating environment has changed.. 

Additionally, process 80 determines 90 whether the value of the optimization function 
for the "better" guess corresponds to a global extremum in the presence of a distortion (i.e., 
whether the value belongs to Sg-d). If the value belongs to set Sg-d, process 80 provides 92 a 
warning to the user and also registers 88 the new guess for the probe's gain, location, and 

25 orientation for the viewer to see. For example, the warning may be an audio signal or a 
flashing signal on a computer display for a user to hear or view. 

Further, process 80 also determines 94 whether the value of the optimization function 
corresponds to a local extremum (i.e., whether the value belongs to S L ). If the value belongs 
to S L , process 80 determines 96 whether the loop count (LC) is greater than a preselected 

30 timeout value (LC max ). If LC > LC^x, process 80 generates 98 a timeout warning. If LC < 
LC max , process 80 loops back 99 to generate a better new guess for the probe's gain, location, 
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and orientation (i.e., process 80 ignores the current new guess). In the loop back 99, process 
80 selects a new initial guess for the probe's gain, location, and orientation (e.g., by randomly 
selecting a new gain (g) and a new point (x, y, z, q>, 8) in the probe's coordinate space). 
Other processes for selecting the new initial guess so that the process 60 subsequently 
5 produces a better guess for a "global" extremum are known to persons of skill in the art. 

In some embodiments, process 80 tries to obtain a better guess for the gain, location, 
and position of the probe after producing a warning that a distortion is present. For example, 
if the number of field measurements is greater than the number of parameters plus one, 
process 80 may produce a better guess by discarding one of the measured field values Bi_ n 

10 measured repeating process 80. If the distortion affects only the discarded field value, 

discarding the distorted value will produce a better estimate of the probe's gain, location, and 
orientation. One field measurement may be distorted because of the presence of a conductor 
near one field sensor or because of a hardware failure in one sensor. 

Some embodiments of process 80 handle overlapping extrema values, which belong 

15 to more than one of the sets of S L , S G -nd, and S G -d» differently. For extrema values in both S L 
and Sg-d, process 80 may provide 92 a warning, indicate 88 the guessed probe location and 
orientation, and try to find a non-overlapping extrema value by repeating steps 82 and 84 for 
a new initial guess of the probe's gain, location, and orientation. For extrema values 
belonging to the set Sl-Slo, the process may generate a warning that identifies the extrema 

20 value as overlapping and then reselect an initial guess for the probe's gain, location, and 
orientation and re-execute process 60 to try to find a value not belonging to Sl-Slo. Of 
course, overlapping subsets such as S l -Slo may also be empty. 

It is important to note that the system gain factor can also be used as a means for 
detecting measurement distortion. The basis for this is the assumption that the system gain 

25 factor should remain essentially constant and stay within a defined range (which is typically 
determined during system calibration). Deviations from this range are indicative of a change 
in the probe's environment or in the system itself. This system gain factor range (SGFR) has 
to be chosen so that it can accommodate normal system gain factor variations for different 
probe positions and orientations throughout the operating volume. This range may also take 

30 into account gain variations caused by manufacturing tolerances and environmental factors 
(e.g., temperature variations, the presence of conductive objects, etc.). 
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If uniform damping of the magnetic field(s) (as sensed by the magnetic field sensors) 
occurs due to the presence of conductive objects within the operating volume, this uniform 
damping is automatically compensated for when the system gain factor is determined. 
Further, this compensation occurs without any loss in the positional and/or orientation 
5 accuracy of the system. 

Accordingly, this automated uniform damping compensation allows for the 
determination of the probe's position / orientation within metallic tubes (e.g., a biopsy needle 
with a field sensor in its stylet). This allows for the tracking of the tip of a biopsy needle as it 
travels through the human body. This is an important advantage, as biopsy needles tend to 

10 bend and/or flex when passing through the body and, accordingly, the tip's position cannot 
be accurately interpolated. These advantages extend to other metallic tube objects, such as 
endoscopes, brachytherapy applicators etc. 

The system can determine whether the probe is inside of a metallic tube by defining 
the SGFR such that system gain factor values for sensors within metallic tubes are not 

15 included in this range.FIG. 6 is a flow chart for a calibration process 100 that finds values of 
the optimization function at extrema to define the membership of sets Sq-nd and Sl- If 
iterative process 60 (of FIG. 4) generates a real global extremum, the probe's gain, location, 
and orientation found by process 60 correlate closely to the actual gain, location, and 
orientation of the probe. If iterative process 60 generates a false or local extremum of the 

20 optimization function, the probe's gain, location, and orientation found by process 60 do not 
correlate closely to the actual gain, location, and orientation of the probe. 

To perform this calibration, process 100 positions 102 the probe at a selected location 
and orientation. Further, an initial gain is selected for the probe. The initial gain is based on 
an accepted range of gain values for a properly functioning probe. During calibration, the 

25 probe mounts on a mechanical positioning frame (not shown) that positions the probe at the 
selected location and orientation. The mechanical positioning frame is made of materials 
that do not distort magnetic fields and provides a separate measurement of the selected actual 
location and orientation of the probe. The separate measurements may be optical 
measurements or mechanical measurements. Process 100 measures 104 the magnetic field 

30 values that correspond to the selected probe location and orientation. 
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Process 100 selects 106 an initial guess of the gain, location and position of the probe 
for optimization process 60. From the measured field values and the initial guess, process 
100 performs 108 iterative optimization procedure 60 to obtain a better-guess value of the 
probe's gain, location and orientation. The optimization procedure 60 also provides a value 
5 for the optimization function that corresponds to the better guess for the probe gain, location 
and orientation and is a value of the optimization function at an extremum. 

Process 100 compares 1 10 the better-guess and actual coordinates of the probe's gain, 
location and orientation to determine whether both coordinates are in close proximity to each 
other. The better guess (x N , Yn, z Nf (p N , 0 N , g) and actual values (x, y, z, q>, 0, g) of the probe 

10 coordinates are close if the values are within a preselected range of each other in a 

component-by-component fashion. If the better-guess and actual coordinates are close, the 
process 100 labels 1 12 the corresponding value of the optimization function as a value 
belonging to Sg-nd- If the better-guess and actual probe coordinates are not close, process 
100 labels 1 14 the corresponding value of the optimization function as a value belonging to 

15 S L . 

To classify values of the optimization function at each extremum, process 100 loops 
back 116 and repeats steps 106-1 14 for other initial guesses of the probe f s gain, location and 
orientation. These repetitions for different initial guesses cover the whole space (x, y, z, cp, 0, 
g) of possible probe coordinates representatively (e.g., by randomly selecting points in the (x, 

20 y, z, q>, 0, g) space). 

Process 100 also repeats the classification of values of the optimization function at 
extrema for other selections of the actual probe gain, location and orientation. The 
repetitions for other actual probe gain, locations and orientations cover a representative 
portion of the parameter space (x, y, z, 9, 0, g) representatively (e.g., with randomly selected 

25 points). The representative portion may be a portion of the entire space (x, y, z, cp, 0, g) that 
is related to other portions of the space by symmetry rotations. 

These repetitions may generate different extrema values of the optimization function 
that belong to Sg-nd and/or S L . For the least-squares sum, the optimization function has 
smaller values at global minima (i.e., in the presence or absence of measurement distortions) 

30 than at local minima. 
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FIG. 7 is a flow chart for a calibration process 120 that finds values of the 
optimization function belonging to Sg-d- The values of S G -d correspond to extrema of the 
optimization function that occur when distortions affect the iteration process 60 (of FIG. 4). 
Process 120 may be performed separately for each type of distortion that can affect either 
5 magnetic field measurements or the processing of magnetic field measurements. The 

distortions may result from nearby conducting or ferrous objects, nearby field sources, sensor 
hardware failures, field source hardware/software failures, and/or software measurement 
processing failures. 

Process 120 physically sets up 122 a selected type of distortion for the system. For 
1 0 example, the setting up of a distortion may include placing a pair of conducting scissors 

inside the monitored volume or causing a hardware failure in the electronics module. After 

setting up the distortion condition, process 120 positions 124 the probe with the mechanical 

positioning frame and receives the value of the probe's actual location and orientation. 

Further, an initial gain is selected for the probe. The initial gain is based on an accepted 
15 range of gain values for a properly functioning probe. The process also measures 126 the 

magnetic field values that depend on the probe's actual gain, location and orientation. 

Process 120 also selects 128 an initial guess for the gain, location and orientation of the 

probe. 

From the measured field values and selected initial guess, process 120 utilizes 130 the 
20 iterative process 60 to obtain a better guess for the probe's gain, location and orientation. 
Iterative process 60 returns an associated extrema value of the optimization function that 
corresponds to each better guess. Process 120 determines 132 whether the new value is 
better than other extrema values of the optimization function, which were generated from 
different initial guesses for the probe's gain, location and position. For the least-squares sum, 
25 the best extrema value is the minimum value. 

If the new value is better than values associated with earlier accepted guesses, the 
process 120 labels 134 the new value as a value of the optimization function at a global 
extremum (i.e., as a member of Sg-d). Values of the optimization function that belong to S G - 
d indicate the presence of distortion. If the new value is not better than earlier values of the 
30 optimization function associated with accepted guesses, process 120 labels 136 the new value 
as a value that corresponds to a false or local extremum. After classifying the extrema value, 
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process 120 loops back 138, 140 to repeat the search for values of the optimization function 
at other extrema by selecting a different initial guess for the probed gain, location and 
orientation. The best extrema value for various initial guesses of the probe's gain, location 
and orientation provides a value in Sg-d. 

5 The process also repeats the search for values of the optimization function in Sg-d for 

different actual gains, locations and orientations of the probe. For each actual gain, 
orientation and location, an extrema value of the optimization function in S G -d may be 
generated. Similarly, different positions of a distorting object (e.g., objects 30 and 32 of FIG. 
1) may generate different extrema values of the optimization function, which belong to Sg-d- 

10 In some embodiments, the extrema values of the optimization function for different 

types of distortion are distinguishable. Distinguishable extrema values fall in different 
ranges. In such embodiments, calibration process 120 is performed separately for the 
different types of distortion to obtain ranges of extrema values of the optimization function 
for each separate type of distortion. Process 80 (of FIG. 5) uses the extrema values to 

15 classify the type of distortion, e.g., hardware failure, software failure, or nearby conducting 
object. 

FIG. 8 shows a computer 142 that determines the gain, location, and orientation of a 
probe from measurements of magnetic fields and indicates the presence of distortions to the 
determinations. Computer 142 may be an embodiment for computer 28 of FIG. 1 or 

20 computer 54 of FIG. 3. 

Computer 142 receives data on measured magnetic field values from a line 144 that 
connects to an output of an electronics module 146. The module 146 may be module 26 
(FIG. 1) or module 52 (FIG. 3). Computer 142 processes the data according to processes 60 
(FIG. 4) and 80 (FIG 5) to determine a probe's gain, location, and orientation and the 

25 presence or absence of distortion. Computer 142 displays the results of the determinations 
on a screen 148. Computer 142 has a data storage drive 154 (e.g., a hard disk or optical 
drive) for storing the above-described executable programs / processes and the data generated 
therefrom. 

A number of embodiments of the invention have been described. Nevertheless, it will 
30 be understood that various modifications may be made without departing from the spirit and 
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scope of the invention. Accordingly, other embodiments are within the scope of the 
following claims. 
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